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We propose a scheme in which broadband nanostructures allow to generate squeezed light and 
entanglement of quantum emitters that are extremely far detuned. It is shown that the reduced 
fluctuations of the electromagnetic field arising from collective resonance fluorescence provide also 
a means to detect the entanglement between the emitters. Due to the near-field enhancement 
in the proposed hybrid systems, these nonclassical effects can be encountered outside both the 
extremely close separations limiting the observation in free space and narrow frequency bands in 
high-Q cavities. Our approach permits to overcome the limitations of noninteracting single emitters 
and is more robust against phase decoherence induced by the environment. 
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The generation of entanglement, one of the most inter¬ 
esting and nonintuitive resources in quantum optics, lies 
at the heart of every physical process that surpasses clas¬ 
sical limits. Entanglement between quantum emitters is 
necessary to perform nonclassical information processing 
fT and to go beyond the standard sensitivity limits in in¬ 
terferometry [2]. Nonclassicality can also be encountered 
in squeezed states of light | 3 ] and becomes manifest in 
the reduced quantum fluctuations of the electromagnetic 
field below the shot-noise level that provide the key to the 
realization of quantum-enhanced applications in commu¬ 
nication, spectroscopy [ 4 ], and imaging [ 5 ]. The close 
connection between different manifestations of quantum 
correlations such as entanglement and squeezed states 
in spin-ensembles [2], also opens interesting questions re¬ 
garding the transfer of entanglement in matter to squeez¬ 
ing in light mm in order to share its nonclassical prop¬ 
erties for quantum applications [8]. Despite extensive 
studies of the creation of entanglement between emitters 
by means of continuous-wave squeezed light [91 HO], the 
possibility to characterize entanglement via the detection 
of squeezed light remains widely unexplored. 

So far, the emission of squeezed light in resonance flu¬ 
orescence mm and the generation of entanglement be¬ 
tween quantum emitters m have been demonstrated in 
a fundamental manner, mainly in free space or micro¬ 
cavities [TO] , Recently, there is a strong interest in ex¬ 
ploring these quantum effects in new regimes in the cou¬ 
pling strength between matter and light by bringing emit¬ 
ters to integrated nanophotonic systems [H[I 5 ]. The 
strength of nanostructures relies on their ability to mod¬ 
ify both the electromagnetic near and far-field [ 16 ], which 
allows for enhancing and controlling the interaction be¬ 
tween emitters 03 - Several schemes for such hybrid en¬ 
vironments have been suggested for the generation of en¬ 
tanglement between quantum emitters m H3 HBH22], 
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whereas only recently a nanoscopic source of squeezed 
light has been proposed [ 23 ] . At present, the experimen¬ 
tal observation of entanglement in nanophotonic systems 
remains elusive. One remaining difficulty is that most of 
the proposed mechanisms rely on either interference aris¬ 
ing in detection with low generation probabilities EU, or 
on the availability of almost identical emitters [T8Tf20] . 
The latter is notoriously difficult in common solid state 
emitters, as local properties of the host matrix inhomo- 
geneously distribute the optical transitions, typically by 
several GHz [ 24 ff 26 ] . These lie outside the reach of nar¬ 
rowband cavities and thus require sophisticated tuning 
methods. 

In this Letter we show how broadband nanostruc¬ 
tures assist in the generation of squeezed light arising 
from the deterministic entanglement of two far detuned 
emitters. We identify the connection of the entangle¬ 
ment between the emitters and the reduced light fluc¬ 
tuations in cooperative resonance fluorescence assisted 
by nanoarchitectures. Our scheme is based on a two- 
photon process in a coupled pair of two-level emitters 
with transition frequencies separated by many linewidths 
(u) 2 —wi = S > 7). Fig. & illustrates the four-level struc¬ 
ture of the coupled system. Two-photon transitions from 
the fundamental state | G) to the doubly excited state 
| E) at the midfrequency (uq + cj2)/2 (blue arrows) be¬ 
come allowed as a result of the dipole-dipole coupling 
potential T1PI12. For emitters coupled by the extreme 
near field in a bulk medium [ 27 ] . incoherent nonlinear 
transitions have been verified experimentally m- Here, 
we exploit the coherent regime of the transition. This 
is possible when PI12 becomes comparable to the emit¬ 
ters’ decay rate 7 and when the driving fields, charac¬ 
terized by the Rabi frequency Pi and the laser detuning 
A = ujl~ (<^i+^ 2)/2, are weak. The coherence, however, 
degrades in stronger fields, when the doubly excited state 
is populated via collective transitions (oc £2 4 ff 12 ) or - to 
a lesser extent - single-photon transitions (oc ff 4 A) [ 27 ] . 
Enhanced dipole-dipole coupling mm is in reach out¬ 
side the extreme near-field regime of bulk environments, 
in state-of-the-art subwavelength broadband cavities [ 29 ], 
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FIG. 1: a) Level scheme of the coupled system, consisting 
of the fundamental state | G) and doubly excited state \E) 1 
connected by a nonlinear two-photon process A = ujl — (w i + 
cj 2 )/2 = 0 (blue arrow). Due to the emitter detuning 5 and 
coupling Qi 2 , the dark and bright single-excitation states \A) 
and |S') are nondegenerate, b) Schematic setup consisting of 
two emitters coupled to a nanoantenna and illuminated by a 
coherent driving field. The observation takes place in the far 
field, c) Incoherent fluorescence signal / = JT(.ET • Ef) as a 
function of the driving field detuning A and free-space Rabi 
frequency Qo- The configuration of two detuned emitters (A = 
780 nm, S — 4007 o,^i 2 — — 6 . 470,712 — — 2 . 670,7 — 2 . 970 ) 
corresponds to distances a — 25 nm from a R — 40 nm gold 
nanosphere. The local driving held enhancement near the 
nanosphere by a factor Q « 2Qq is taken into account. 


nanowaveguides [14, [25|, or antenna structures 15 . All 
these hybrid systems can provide the necessary interac¬ 
tion to overcome saturation broadening and permit much 
larger emitter separations or even far-held coupling. We 
hnd that the nanostructure-assisted nonlinearity allows 
for the generation of entanglement and, in consequence, 
the emission of squeezed light, exceeding the performance 
of uncoupled emitters. 

As a proof of principle, we consider the system de¬ 
picted in Fig.[l] 3 , consisting of two quantum emitters 
placed at a distance a = 25 nm from a sub wavelength 
gold nanosphere of radius R = 40 nm and symmetrically 
driven. Each of the emitters (i = 1, 2) is characterized by 
the coherence operator of a two-level emitter, = \g)(e\ 
with resonance frequencies separated by S = 400 70 . Near 
the nanostructure, not only the dipole-dipole coupling 
f^i 2 and the incoherent coupling rate 712 , but also the 
emitter linewidth 7 and Rabi frequency 12 are altered 
from their values in free space, 70 and 12 o, increasing the 
bandwidth of the interaction and requiring lower driving 

powers. The positive frequency helds E i = gfti, scat¬ 
tered by each of the emitters in the presence of the nanos- 
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FIG. 2: Steady state concurrence as a function of the driving 
field amplitude and detuning. Parameters as in Fig. [l] Bipar¬ 
tite entanglement is excluded in the white area for detuned 
emitters (A > Q 12 ) and is restricted to the region below the 
white dashed curve for identical emitters (A <C ^ 12 ). 


tructure to the far-field, can be expressed in terms of the 
classical Green’s tensor. We then obtain the steady-state 
correlations of the total scattered field from the quantum 
optical master equation (see Supplemental Material for 
details). 

Since the coherent dynamics are limited by the pop¬ 
ulation of the excited states, we first characterize the 
response to the driving field in terms of the incoherent 

fluorescence signal I = J2 i (E i -E i ), which is a common 
experimental measure of the excitation spectrum and al¬ 
lows to identify the dipole-dipole coupling [24] . For rela¬ 
tively equivalent positions of the emitters and a balanced 
detector position, we may assume both equal scattering 
amplitudes and decay rates 7 ^ = 7 (the impact of asym¬ 
metric configurations is discussed in the Supplemental 
Material). Fig.JlJ displays the incoherent fluorescence as 
a function of the laser detuning A and the free-space 
driving amplitude 12 q. Its onset suggests a rough bound¬ 
ary of the coherent regi me. Apart fr om the single-photon 
resonances at A « ± yf(S/2) 2 + there is a clear non¬ 

linear absorption peak at A = 0 that arises from the 
direct two-photon transition from the fundamental state 
| G) to the doubly excited state | E). The peak saturates 
for stronger driving 12 , where it eventually overlaps with 
the saturation-broadened single-emitter lines. The collec¬ 
tive fluorescence peak at A = 0 would not be observable 
at comparable emitter separations and detuning in free 
space since the required stronger excitation field would 
mainly increase the single-emitter background. More¬ 
over, the intensity observed in the far field (oc \g\ 2 ) is en¬ 
hanced by a factor ^1.7 with respect to free space [23] 
and it can still be significantly improved in engineered 
structures [HI 3jJ. We finally notice that the large emit¬ 
ter detuning considered here would usually lie outside the 
interaction range of solid-state emitters in ultrahigh-Q 
cavities. For example, for cryogenic molecules with typ¬ 
ical optical linewidths 7 o/27t ^ 50MHz [24] . the equiva¬ 
lent detuning S/2tt = 20 GHz would barely lie within the 
bandwidth of a cavity of Q = 10 4 . 

The coherent excitation via the nanostructure-assisted 
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dipole-dipole interaction makes it possible to generate 
steady-state entanglement between the two emitters de¬ 
spite their large detuning. We identify bipartite entan¬ 
glement in the following by nonzero values of the concur¬ 
rence (7, which ranges between 0 for fully separable states 
and 1 for maximally entangled states ED- Fig-0 shows 
an evaluation of the steady-state concurrence as a func¬ 
tion of the driving field detuning and amplitude for the 
previous configuration. A significant amount of entangle¬ 
ment can be reached close to the two-photon resonance. 
As Q is increased, the onset of saturation in the collec¬ 
tive transition shifts this maximum in frequency (asym¬ 
metrically due to unequal populations of the entangled 
states | A) and |S')), and finally suppresses entanglement. 
This behavior differs strongly from the case of identical 
emitters ((5 <C £2,7), for which single-emitter saturation 
suppresses the generation of entanglement at comparable 
driving £2 (in the region above the dashed white curve in 
Fig©. Moreover, the maximum concurrence C ~ 0.3 is 
comparable to typical values for identical emitters cou¬ 
pled via nanostructures [18]. 

Quantum correlations between the emitters EnnBiig 
strongly affect the statistical properties of the scat¬ 
tered light. These can enhance squeezed resonance 
fluorescence [33] that is characterized by quadrature 
fluctuations below the shot-noise limit and can be ex¬ 
perimentally addressed in homodyne or heterodyne se¬ 
tups 13 H3. For a general electric field quadrature 

E e = e l6 Ei + h.c., we identify reduced fluctuations 
of a single vector component by a negative value of the 
normally ordered variance (AS) 2 = (: (Eq — (Eg)) 2 :) 
= (A£i) 2 + (Af 2 ) 2 + (A£i 2 ) 2 . Here, the fluctuations 
from single emitters correspond to 


( A £;) 2 
21 | 2 



e 2i <»+**) (*.) 2 


( 1 ) 


with scattering phases <pi = arg(^) and recover previous 
results [23 in the absence of interemitter coupling. The 
interaction modifies each of these terms and gives rise to 
the additional cross-correlations 

e i(<Pi-4>2) _ (<t!)(<ti)^ (2) 

+e i(20+0i+0 2 ) (^CTi) - (<7 2 )(<7 i)^ 


(A£ 12 ) 

2|siS2| 


= 2Re 


that clearly vanish for uncorrelated emitters. Therefore, 
negative values arising from the cross-terms denote re¬ 
duced fluctuations that allow to identify bipartite entan¬ 
glement. 

The black curve in Fig.[3^i shows the optimum degree 
of squeezing (AS) 2 / 2 |g| 2 , obtained by varying the phases 
as a function of the laser detuning. The driving 
£2o = 6 O 70 was chosen well below the two-photon satura¬ 
tion. As before, we focus on a favorable balanced detec¬ 
tion configuration (\gi\ = \g\, see Supplemental Material 
for a discussion of the general case). A distinct mini¬ 
mum of the fluctuations arises close to the two-photon 
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FIG. 3: a) Electric field fluctuations (black curve) obtained 
by optimizing the detection quadrature angle and steady state 
concurrence (red curve) vs. laser detuning at £2o = 6 O 70 . The 
gray shading indicates squeezed light. The green and blue 
dashed curve gives the effective two-level approximation on 
the | E) — | G )-transition for the squeezing, and the concur¬ 
rence, respectively, b) Electric field fluctuations as a function 
of the driving field detuning and amplitude, optimizing the 
detection quadrature angle. No squeezing can be observed 
above the black dashed limit. The red curves indicate the 
limits for independent emitters. The fluctuations exceeds the 
minimum value feasible from two independent emitters inside 
the white contour. Parameters are chosen as in Fig. [l] 


resonance A = 0 and mirrors the behavior of the concur¬ 
rence, indicated by the red (upper) curve. To exemplify 
the link between entanglement in matter and squeezing 
in light and to understand the main mechanism under¬ 
lying the present scheme, we consider a simplified model 
for weak driving on the two-photon resonance. An ef¬ 
fective two-level picture arises, because the single-photon 
transitions are far detuned and the dipole-dipole coupling 
mainly induces a coherence peg between the fundamen¬ 
tal and doubly excited state, so that the density matrix 
is dominated by this element and the populations. From 
these elements, we find from Eqs. Q and ([2| the normal¬ 
ized degree of squeezing 


(A £) 2 

2|«?| 2 


1 — Pgg +Pee 


e ^+^+^]p EG 


(3) 


which is optimized when the last term equals —\peg\- 
The approximation (green dash-dotted curve) agrees ex¬ 
cellently with the full squeezing amplitude (black curve). 
Interestingly, Eq. © is closely related to the spin squeez¬ 
ing (2113 present in the system and reinforces the link 
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FIG. 4: Maximum degree of squeezing as a function of the 
driving field detuning and amplitude, obtained by optimizing 
the detection quadrature angle in the presence of strong ad¬ 
ditional pure dephasing processes at a rate 7 * = 27 for which 
single-emitter squeezing is fully suppressed (cf. red dashed 
limit). All other parameters as panel Fig. [I] Inset: Analo¬ 
gous representation for the concurrence. 


to entanglement. In fact, a similar treatment of the con¬ 
currence near the resonance (blue dashed line) is equally 
dominated by \peg\ ( see Supplemental Material for de¬ 
tails). Altogether this identifies the two-photon coher¬ 
ence as the source of both squeezing and entanglement 
and supports the description in terms of an effective two- 
level model. 

A full evaluation of the electric field fluctuations as a 
function of the driving field amplitude and detuning is 
given in Fig-Hp At the chosen parameters, the strongest 
squeezing is encountered at Do ~ 75yo near the center 
frequency A « 0 and closely follows the regime that 
allows for maximum bipartite entanglement in Fig.|2j 
When the driving field intensity is further increased, sat¬ 
uration of the two-photon transition suppresses squeezing 
at A = 0, similar to the concurrence. As the collective in¬ 
teraction also affects Eq. 0 , the boundaries of squeezed 
light can surpass the ones for uncoupled emitters [23] 
(denoted by the red dashed limits in Figj3k). The 
strongest squeezing with values (AS) 2 / (2\g\ 2r ) « —0.21 
significantly exceeds the universal threshold —0.125 that 
limits the squeezed resonance fluorescence from indepen¬ 
dent emitters m and is indicated by the white contour 
in Fig. [ 3 ) 3 . Thus, a strong collective enhancement of the 
squeezing per emitter [33] can be achieved despite their 
large detuning. This is a general result for two-level sys¬ 
tems coupled by means of strong dipole-dipole interac¬ 
tions and can be used to identify entanglement. 

As an advantage for the realization of nanoscopic 
sources of squeezed light, the collective effects induced 
by the nanoarchitecture allow for overcoming significant 
additional pure dephasing [34]. This phenomenon oc¬ 
curs commonly in solid-state environments due to the 
coupling to thermal phonons and it presents a main lim¬ 
itation to the operation time of quantum gates. For in¬ 
dependent emitters, additional pure dephasing at a rate 
7 * > 7/2 will fully suppress the generation of squeezed 
light, even though the enhancing effect of a nanostruc¬ 


ture may modify the boundary significantly [23]. At the 
onset of this unfavorable scenario, we represent in Fig. [4] 
the minimized electric field fluctuations as a function of 
the driving field and detuning. The full suppression of 
squeezing for uncoupled emitters is indicated by the red 
dashed limit, i.e. (AS) 2 > 0 for all Do- Still, a significant 
amount of squeezing survives near the two-photon reso¬ 
nance, though at a reduced absolute value. The optimum 
follows even more closely the behavior of the concurrence 
(shown in the inset) than in the absence of pure dephas¬ 
ing because terms stemming from the single emitters are 
now suppressed. These facts reflect that squeezed emis¬ 
sion identifies entanglement despite additional dephas¬ 
ing and underline that nanostructures allow to taylor 
collective interactions, providing more robust sources of 
squeezed light in the presence of decoherence. 

The combination of nanophotonic structures with 
quantum emitters presents a promising prospect for novel 
integrated sources of squeezed light and the generation of 
multipartite entanglement at the nanoscale. The strong 
local field enhancement in nanoarchitectures feasible over 
a broad frequency range can provide significant intermit- 
ter coupling outside the regimes of extreme near fields in 
bulk media. Using a metallic nanosphere as a paradig¬ 
matic example, we have found that both entanglement 
and squeezing can be achieved for two emitters that are 
mutually far-detuned, beyond the limits that are spec¬ 
trally in reach in typical narrow-band cavities. This can 
greatly relax the necessity of tuning precisely the opti¬ 
cal resonances of solid-state emitters required to achieve 
large degrees of entanglement in hybrid systems and it 
can assist the future detection of entanglement. More¬ 
over, we have shown that the generation of squeezed light 
from interacting quantum emitters can serve as a reliable 
witness of entanglement. Therefore, resonance fluores¬ 
cence from interacting emitters can help to explore the 
transfer of nonclassicality from matter to light and the 
distribution of useful quantum states in photonic devices 
that surpass the limits of classical optics im Our re¬ 
sults show that collective interactions mediated by nanos¬ 
tructures can allow to overcome a significant amount of 
additional pure dephasing, even beyond the improvement 
achievable in independent-emitter squeezing in a nanos¬ 
tructure [23] . As both effects counteract phase decoher¬ 
ence - a fundamental limit in setups at nonzero tempera¬ 
tures - nanostructures can therefore allow to scale up the 
operating temperatures of solid-state emitters. A simple 
estimation for typical cryogenic molecules [35] shows that 
an increase of 300 in the radiative linewidth could allow 
to work at temperatures « 70 K rather than below 5K, 
which would have significant technological implications. 
Due to the generality of the formalism, our results may 
be transferred to other broadband environments such as 
optical antennas, subwavelength cavities, or nanowaveg¬ 
uides, where high coupling strengths may allow to push 
the limits for the generation of squeezed light even fur¬ 
ther. In addition, we expect that the collective effects 
at work can be scaled up to higher numbers of emitters 
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nanni Eg, a future study of which will also provide a 
better understanding of the origin of squeezing in struc¬ 
tured nonlinear media. 
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given by (see, e.g., Refs. Ham]) 

p= ~^[{H - ]T Q, m n^ m ^n),p\ + C[p] , (7) 

m^n 

C iP\ = E ( a m a nP + paln^n ~ ^mP^V) • ( 8 ) 

m,n 

Photon-mediated interactions result in modified dipole- 
dipole potentials Cl mn while the dissipative dynamics, i.e. 
the individual and cooperative decay of dipole excita¬ 
tions, are described by the Liouvillian C. Here, the rates 
are the real and imaginary parts of 

^ran T 17ran/2 — d m ' 9n (**m) • (9) 


A. Master equation in macroscopic 
electrodynamics 

We use the approach of macroscopic quantum electro¬ 
dynamics [37] to describe the behavior of N two-level 
quantum emitters in an arbitrary environment. Each 
emitter is described by its transition dipole element d n , 
emitter operator <r n , and transition frequency uj n . Within 
the dipole, rotating wave, and Markov approximation the 
full positive frequency field operator is given by 


Due to the local field enhancement, the single-body decay 
rates in the nanostructure = 7 mm generally differ 
from the free-space value 70 . For this work, the master 
equation 0 is solved numerically for the steady state 
solution in the case of N = 2 . Aditional pure dephasing 
at a rate 7 ^ can be included phenomenologically [34] by 
adding the term 

^d[p] = ^ ^ / ymi ( ^z,m < ^z,mP T P&z,m&z,m 2<J z ,m P&z ,m) • 

m 

( 10 ) 


N 


£ to t(r, t) = £free(r, *) + £ K ( r > *) ( 4 ) 


n= 1 


with E n (r, t) =g n ( r )a„(t) 


(5) 


Here, g n {r) « ^ G(r,r n ,co 0 ) ■ d n encodes the pho- 
ton propagation in an arbitrary broadband environment 
via the classical dyadic Green’s tensor, where r denotes 
the detector position. As long as dispersion in a nanos¬ 
tructure is weak, we can evaluate the Green’s tensor at 
the center frequency uj 0 = Since we con¬ 

sider far-held observation, we neglect a small position- 
dependent quantum correction connected with dispersion 
potentials, which would require some care in the study 
of near-helds [22]. For a gold nanosphere, we evaluate 
the Green’s function semi-analytically |38l 139] . using a 
Drude-Lorentz interpolation of tabulated optical data for 
gold [401 . 

Under these conditions, the evolution equation is ex¬ 
pressed in the interaction picture with respect to the 
free quantum fields and in the rotating frame of a clas¬ 
sical driving held EL(r)e~ luJLt + h.c. that translates to 
single-emitter Rabi frequencies Vt n = Ei,{r n ) • d n /Ti . The 
Hamiltonian reads 


h = %y j 


K - U L )(T z ,n ~ ( + h.C. 


(6) 


with <j z , n the population inversion of particle m. The 
quantum optical master equation in the Lindblad form is 


To obtain the solution and the relevant correlation func¬ 
tions, we employ the single-body basis | 1 ) = \g1g2), | 2 ) = 
\91e2), |3) = |ei g 2 ), 14) = |eie 2 ), where e* () denotes the 
i-th emitter in the excited (ground) state. 

To describe the physics involving the weakly driven 
two-photon transition, it is insightful to use the collec¬ 
tive coupled eigenstates in the absence of driving. For 
two coupled emitters detuned by S = UJ 2 — ^1 and in 
the absence of driving, the coherent terms in the master 
equation 0 are diagonalized by the states 


\G) = 11), 

1 E) = |4), 

(ii) 

\S) = a 3) + 6 2 ), 

|7L> = 6|2) - o|3) , 

(12) 

d 

n — 

h- n±2 

(13) 

Vd 2 + n 2 u 

Vd 2 + n 2 u 

d = 25+ Ja5 2 + n ; 2 . 


(14) 


The corresponding eigenfrequencies lje = = (wi + 

6 ^ 2 )/ 2 , c os = —ua = v^?2 + 4£ 2 provide the level 
scheme shown in Fig. la in the main text and describe 
well the spectrum encountered at sufficiently weak driv¬ 
ings. In the case of identical particles S = 0,71 = 72 , the 
coupled basis recovers the bright and dark maximally 
entangled Dicke states (\S), \A) |+), |—)) with decay 
rates = 7 d= I 712 I, respectively. These states govern 
the dynamics at weak driving Hang, whereas in the de¬ 
tuned case S 7 ^ 0 , | G) and | E) are mainly responsible for 
the response at the two-photon resonance as discussed in 
the main text. 
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B. Incoherent fluorescence 

One relevant experimental observable used to identify 
the two-photon resonance is the incoherent fluorescence 
as a function of the driving frequency [24], which can be 
measured outside the spectral range of the exciting laser. 
By assessing where the fluorescence level is still low, one 
can identify the parameters for which coherent effects 
are expected to be significant. The fluorescence signal is 

proportional to I = J2i(^i ' &t) = Hi \9i(r)\ 2 (a\<Ji). 
In the main text we set |^i (r*) | = |g 2 (T)|, which is a 
reasonable assumption when the emitters are placed at 
equivalent positions with respect to the nanostructure. 
This can be inferred from Eqs. , where for similar 

emitter positions the one-point Green’s tensor provides 
Gr(ri, ri, cjo) ~ G(r 2 , r* 2 , cjo) and at equivalent far-held 
detection distances G(r, rq, cjo) ~ G(r,r 2 ,^ 0 )- Nev¬ 
ertheless, we have numerically checked that values dif¬ 
fering by up to an order of magnitude between |gi(r)| 2 
and 1 02 (t)| 2 and between 71 and 72 , the impact of the 
two-photon transition remains qualitatively similar to the 
case shown in Fig. 1 in the main text. Mainly, this in¬ 
troduces an asymmetry in the fluorescence and a gives a 
stronger weight in intensity to the emitter with the larger 
scattering amplitude \gi\. On the other hand, the states 
| G) and | E) are still connected by a collective transition 
at A = ujl ~ (cji T CJ 2 )/2 = 0. This can be seen from the 
perturbative formula for the population of the excited 
state deduced from Ref. m 

P44= (7i+72) 2 + (2A)2 X (15) 

^ 2 f 4 i^ 2 (kh + <^2 — ^12[^1 + ^2] 

(cji -u l )(uj 2 -Ml) (cji -uj l )(uj 2 -Ul) 

which can still provide comparable values for 71 7 ^ 72 or 
asymmetric driving conditions 7 ^ Qr ,2 • 


C. Concurrence near the two-photon resonance 

The amount of entanglement generated in the coupled 
system may be quantified for a general quantum state 


using Wootters’s concurrence m 6 = max[{ 0 , Ai — A 2 — 
A 3 — A4}], where the A 2 are the eigenvalues of the operator 
p(<J y , 1 (8) o'y,2)p*(o'y,i ® a yp) in descending order. The * 
denotes the complex conjugate. To gain a better under¬ 
standing of the entanglement induced by the two-photon 
resonance, we consider the limit of a weak driving field. 
In this case, the solution will be close to that of an un¬ 
driven system, for which exact analytical solutions are 
available. Here, the density matrix keeps a cross-shaped 
symmetry |42l and the concurrence is expressed as 

6 ^max{ 0 , 61 , 62 } with 

61 =2|p4i| — 2y/p22p33 ? C 2 = 2|p 2 31 — %y/PllpAA • 

(16) 

We have checked for our scheme that 61 applies near the 
two-photon resonance at weak driving, where the single¬ 
photon populations are negligible. This shows that en¬ 
tanglement is due to the driving of the two-photon co¬ 
herence P 14 connecting the double ground and double 
excited states. When |7 12 1 > 0, the asymmetry between 
the dark and bright states is visible in a dispersive con¬ 
tribution to the concurrence at larger positive detunings, 
for which then 62 dominates. 


D. Electric field quadrature fluctuations 


We identify sub-shot-noise fluctuations in an electric 

field quadrature component Eq = e l ° E + e~ l0 E by 
negative values of its normally ordered variance [4] 

(:(AE e (r,t)] 2 :) = 2Re [e i2e ((E + (r,t)} 2 ) - e™ (E + (r,t)) 2 } 
+ 2(E~(r,t)E + (v,t))-2\(E~(v,t))\ 2 , (17) 

where E = corresponds to the negative- 

frequency electric field. Eq. (|4) links the field variance 
to the emitter correlation functions. We find the elec¬ 
tric field fluctuations for one vector component (index 
suppressed for clarity) 


(• [AEfl(r, t )] 2 :) 

(A£i ) 2 

2M 2 

(A £ 2 ) 2 

2 M 2 

(A£i 2 ) 2 

2151521 


(AS,) 2 + (A S 2 ) 2 + (A £ 12 ) 2 , (18) 

533 + 544 — |5i3 + 5241 2 — Re [e 2l ( e+ ^7p3i + 54 2 ) 2 ] , (19) 

522 + 544 — |5i2 + 5341 2 — Re [e 2l ^ 6+ ^ 2 \pi2 + 534) 2 ] ! (20) 

2Re [e^-^piz + e~^ 29+ ^p u ] ~ 4Re [e~W+^ >(p 13 + p 24 )]Re [e- i(e+M (p 12 + 53 4 )] • 


( 21 ) 
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In the present scheme, two-photon squeezing can be ob¬ 
served best in a balanced detection configuration, where 
\gi\ 2 /\g 2 \ 2 « 1. In very unbalanced situations, the signal 
is dominated by the single emitter terms [Eqs. (19) or 
( |20| )] and reduce the visibility of the cooperative squeez¬ 
ing term [Eq. ©]■ Nevertheless, we have confirmed 
numerically that the cooperative term can be discerned, 
as in the case of the incoherent fluorescence, for differ¬ 
ences between the emitter’s coupling constants up to an 
order of magnitude. 

At weak driving near the two-photon resonance, we can 
derive a simple semianalytical formula following closely 
the approximations provided in the previous subsection 
for the concurrence. We assume again a cross-shaped 
density matrix, for which Eq. (18) reduces to 


-(:[AE e (r,t)] 2 :) « \gi\ 2 (p33 + Paa) + |^ 2 1 2 ( P22 + Paa) 

+ 2\ gi g 2 \Re [e^~^ p 23 + e~ i( - 2e+ ^ + ^p 14 ]. (22) 


Since the single-photon coherences are small near the 
two-photon resonance at weak drivings (|p 23 1 |pi 4 1), 
we can further simplify the squeezing signal. Thus for a 
balanced situation \gi \ = \g 2 \ = |#|, 0 i = 02 we obtain 


where we have redefined the quadrature angle 0' = 0 + 
(01 + 0 2 )/2 without loss of generality. Reduced electric 
field fluctuations are thus traced back to the two-photon 
coherence pAi and are suppressed by any depopulation of 
the ground state and optimized for a quadrature angle 
such that —2Re [e~ 2ie ' pia] = — 2|p 14 |. 

It is instructive to point out the connection to spin¬ 
squeezing, considering the fluctuations of the collective 
operators S x = E„(<4 + <?n)/2, S y = E„(<4 ~ n)/2i, 
and S z = ^2 n cr z ^ n . For driving near the two-photon 
resonance, we find that (S x ) « (S y ) ~ 0. Spin squeezing 
M occurs when 


2((5g) ~ (^) 2 ) 

1(^)1 


< 1 . 


(24) 


Direct algebraic manipulation 
this is actually equivalent to 
(:[A5] 2 :} = ((: S 2 :) — (S x ) 2 ) < 0. Here, 


shows that 
the condition 


4( : [AS] 2 :) — 1 — pn + Paa ~ 2Re (pia) 

+ [P 22 + P 33 - 2Re (P 23 )] • (25) 


(:[AE g (r,t)] 2 :) 

2|^| 2 


« P44 + (1 - P 11 ) ~ 2Re [e 2l6 ' p 14 \ 

> -2\p u \ , (23) 


Since the antisymmetric state is not driven, the last 
paranthesis vanishes recovering an expression that is pro¬ 
portional to Eq. (23) for vanishing quadrature phase and 


equal values of gi = g 2 . 
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